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SUMMARY 


A study is made of the steady laminar flow of a compressible viscous 
fluid in a circular pipe when the fluid is accelerated by an axial body 
force . The application of the theory to the magnetofluidmechanics of an 
electrically conducting gas accelerated by electric and magnetic fields 
is discussed. Constant viscosity, thermal conductivity, and electrical 
conductivity are assumed. Fully developed flow velocity and temperature 
profiles are shown, and detailed results of the accelerating flow develop- 
ment, including velocity and pressure as functions of distance, are given 
for the case where the axial body force is constant and for the case where 
it is a linear function of velocity. From these results are determined 
the pipe entry length and the pressure difference required. 


INTRODUCTION 


An electrically conducting fluid, such as a high-temperature ionized 
gas, or a "seeded" gas, flowing in a channel or pipe can be accelerated 
by an electric body force which is the result of application of mutually 
perpendicular electric and magnetic fields at right angles to the direction 
of flow. A study of the electromagnetic interaction with a plasma is of 
interest because of the possibility that the velocity of gas flow in a 
wind tunnel or in the exhaust of a rocket or a jet engine can be substan- 
tially increased by such an electromagnetically induced force . This 
concept is mentioned, for example, by von Karman (ref. l) , and devices 
which utilize it have been investigated and described by Patrick (ref. 2 ) 
and Ghai (ref. 3 ) among others. Thrust augmentation by this method may 
be particularly useful in interplanetary flight because the achievable 
specific impulses are likely to be high. 

The theoretical magnetofluidmechanic channel flow has been treated 
by numerous investigators. In particular, viscous incompressible flow 
and accelerating inviscid compressible flow have been studied extensively. 
Resler and Sears (refs. 4 and 5) have discussed and analyzed in some 
detail certain aspects of fully developed viscous incompressible flow 
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and accelerating inviscid compressible flow. The problem of determining 
entry length for the incompressible viscous flow of liquid metals in 
pipes entering a region of uniform transverse magnetic field has been 
studied by Shercliff (refs. 6, 7 , and 8). Shercliff solved an unsteady 
flow problem for a "settling time" or time to reach the steady state 
after a magnetic field is applied. He then converted the settling time 
to an approximate steady-flow entry length by multiplying the settling 
time by the mean velocity of flow in the pipe . He showed by experiments 
that the correct order of magnitude for entry length had been estimated. 
References 6 and 7 dealt with flow between infinite parallel planes and 
flow in rectangular pipes in which only nonconducting walls were consid- 
ered. In reference 8 flow in a circular pipe with conducting walls was 
studied. For that case the problem was treated in two stages: In the 

first stage thin boundary layers developed with a core of uniform velocity; 
then, in the second stage, the core profile was changed by electromagnetic 
forces to its final form, with the boundary -layer thickness remaining 
unchanged. The field of magnetohydrodynamics has recently been reviewed 
by Sears (ref. 9 ) and "by Rossow (ref. 10) . 

The present study treats the steady laminar flow of a viscous fluid 
accelerated by an axial body force in a circular pipe but differs from 
the study by Shercliff (ref. 8) in that the flow is considered to be 
compressible rather than incompressible . The fluid is considered to 
have a uniform velocity profile as it enters the pipe and to accelerate 
until it approaches the "fully developed flow" condition, in which the 
flow is exactly parallel and where there is no further acceleration or 
change in the velocity profile. (For a liquid flowing in a pipe the fully 
developed laminar flow is the familiar Hagen-Poiseuille pipe flow, 
discussed in reference 11. Some effects of compressibility, variable 
properties , and body forces on fully developed laminar flows in two- 
dimensional channels have been studied by Maslen (ref. 12).) It is shown 
that, when fully developed flow of a compressible fluid is defined simply 
as a nonaccelerating parallel flow (which definition applies also to incom- 
pressible fully developed flow), the necessary conditions for such a flow 
to exist are that the axial gradients of density, pressure, and temperature 
be zero and that a body force be supplied to sustain the flow. When these 
conditions are supplied, via proper heat transfer to the wall and, for 
example, electromagnetic body forces, then the body forces and viscous 
shear forces can become balanced and the flow will approach an asymptotic 
state . 

Solutions for the fully developed flow velocity and temperature 
profiles are found both for the case where the axial body force is 
constant and for the case where it is a linearly decreasing function of 
velocity. The development of the velocity profiles in the accelerating 
portion of the flow is treated approximately. It is shown how these 
solutions can be applied to the study of the magnetogasdynamic pipe -flow 
problem in which the fluid is a conductor of electricity and the body 
force is induced by electric and magnetic fields. 
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SYMBOLS 

a 

radius of pipe 

A 

defined by equation ( 15 ) 

An 

Fourier coefficient 

b i 

Jy H 0 

0 in case I (see eq. ( 9 )) 

(TUqHo 2 

A 

3 B n 

9 11 

% 

■ ■ in case II 
H 0 -u 0 

6 

B n 

Fourier coefficient defined by equation (A3) 

Cl, c 2 

constants of integrations 

C 

defined by equation (^6) 

C P 

specific heat at constant pressure 

B n 

Fourier coefficient defined by equation (A6) 

E 

electric field 

* 

E e 

electric field due to external effects 

E ± 

electric field due to internal separation of charges 

0, 

Joule heating energy (eq. ( 5 ^)) 

a 

8,2 n 

PU 0 2 

f 

function of A n; defined by equation (AlO) 

F 

axial body force per unit vc^Lume 

— > 
F 

body force vector per -unit volume 

F 

a 2 

g 

function of A n , defined by equation (All) 

G 

function of 3c, used first in equation (4o) 

h 

function of A n , defined by equation (A17) 

-► 

H 

magnetic field vector 
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Ho 

Io, lx 
J 

Jo J Jl 
k 

K 


applied magnetic field in the z direction 

modified Bessel functions of the first kind of orders 
zero and one 

current density 

Bessel functions of the first kind of orders zero and one 
thermal conductivity 


g H 0 2 a 2 


, Hartmann number 


Mo 

m 

m 

n 

No 


inlet Mach number 


pu 


- m = N 0 pu 

index of eigenfunctions 


Po u o a 


, Reynolds number at pipe inlet, based on radius 


P 

P 

Pr 

— > 
q. 

r 

r 

R 

Rn 

Si, S2, S3 

T 

T 


thermodynamic pressure 
a 

P^o P 


PC. 


•; Prandtl number 


velocity vector 

distance normal to pipe axis 
r 

— m 

a 

gas constant 

function of r; solution of equation ( 28 ) 

defined by equations (Aik) , (A15) , and (A20 ) , respectively 

temperature 

CpT 

Un 2 
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u 

u 

U 


v v 


X 


X 


X, 


n 


y 


z 

a 

3 

7 

6 n 

'H 

P 

£ 

I 

P 

P 


velocity in the x direction 

u 

u o 

difference in dimensionless velocities, defined by 
equation (2 0) 

dimensionless specific volume, 4 

P 

radial velocity component 

axial distance from the pipe inlet 

x 

a 

function of x 

vertical coordinate perpendicular to pipe axis 
horizontal coordinate perpendicular to pipe axis 
constant of integration (eq. ( 29 )) 
constant in equation (A7) 
ratio of specific heats 

defined by equation (27) for case I; by ( 51 ) for case II 
function of r defined by either equation ( 67 ) or (82) 
coordinate angle in the y-z plane 
eigenvalue 
fluid viscosity 

function of r, defined by equation ( 36 ) 

function of r, defined by equation (A 5 ) 

mass density of fluid 
P 

Po 

electrical conductivity 
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Subscripts and Superscripts 


( 

( 

( 

( 

( 

( 

( 

( 

( 


) 

) 

)’ 

)* 


value at x = o, r = o 
value at large x 
case I (jy = constant) 
case II (Ey = constant) 
value at the wall (r = l) 
component in the y direction 

dimensionless quantity, except where defined otherwise 
vector 

differentiation with respect to r 

value when velocity has approached to within 1 percent of 
fully developed flow velocity on the center line 

EQUATIONS OF MOTION 
Definition of the Problem 


The steady high-speed compressible flow of a continuum gas in a 
semi-infinite straight circular pipe is considered. The fluid enters 
the pipe with high velocity at x = 0 (see fig. l) and each fluid element 
is thereafter acted on by an axial body force . The Reynolds number 
pua/q is assumed to be low enough that viscous forces are important and 
that the flow is laminar. A compressible fluid flowing through a pipe 
may or may not approach an asymptotic state . If conditions are such that 
the viscous shear, pressure forces, and body forces in the fluid can 
become balanced, then the flow approaches the fully developed condition. 
Only cases in which the flow does approach this condition will be studied 
here. The necessary conditions for fully developed flow to occur (as 
described in the Introduction) will be derived in a later section. 


Hypotheses Used and Their Effects on the Governing Equations 


In order to conveniently study the flow, the following hypotheses 
are made: (l) the r and Q components of velocity (fig. l) are nearly 
zero; (2) the flow is axially symmetric j (3) only axial body forces act 
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on the fluid; and (^) the viscosity is constant. The equations which 
result from application of these hypotheses to the flow conservation 
equations must he treated as approximations. 

The first hypothesis indicates that the flow is approximately 
parallel. Although this is exactly true only for the fully developed 
flow, it is taken as a reasonable assumption for the accelerating portion 
of a high velocity compressible gas flow in a straight pipe. The conti- 
nuity equation for the compressible flow of gas is then 

1- (pu) ~ 0 (l) 

ox 

It must be noted that this equation must be treated as an approximation 
to the actual continuity equation; only in the special case of nonaccel- 
erating flow can it exactly describe the true physical flow. For equation 
(l) to be exactly true in the accelerating flow, very special conditions 
on the heat transfer through the fluid and to the walls would be required. 
However, as will be seen, the analysis is greatly simplified by use of 
equation (l), and indications of the accuracy of the results will be 
shown later. The approximate equation (l) shows that pu can be a 
function only of r, not of x; that is. 


pu ~ m(r) 


( 2 ) 


The implication is that each cylindrical lamina of fluid, at a distance r 
from the center of the pipe, moves as in a one -dimensional flow with pu 
constant, the motion in each lamina being affected by the viscous forces 
due to the shear stresses on the adjacent fluid laminae. Hear the pipe 
axis the viscous forces are small and the velocity of the fluid particles 
is greatly increased by the body forces, whereas near the wall, the higher 
viscous forces have a greater retarding effect. The density must therefore 
be greater near the wall. This is analogous to a compressible boundary 
layer adjacent to a cold wall, in which the density is higher near the 
wall. It might be noted that, in compressible boundary-layer solutions, 
one must consider v r because the flow downstream of a point depends on 
what occurred at the point, that is, on the magnitudes of both components 
of velocity. If v r were neglected in a boundary-layer problem where 
the downstream flow is not known but depends on the flow development 
upstream, errors in the solution of the downstream flow would be compounded 
because of inexactness of the solution upstream from a point. However, 
in the present problem, precise knowledge of v r is not necessary to 
prevent such a compounding of errors because the solution is forced to 
approach the known exact solution of the equations governing the downstream 
fully developed flow, as will be shown later. 
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The momentum equations which govern the motion 
compressible Navier-Stokes equations in cylindrical 
following form: 

8u _ 8p 4 8 2 u M- 8 f 
PU 8x " " + 3 ^ Sx 2 + r 8r V 


reduce from the 
coordinates to the 



(3) 


8p 1 8 2 u 

— =- -f- — |JL - 

dr 3 dx dr 




where F is the axial body force per unit volume . If F is constant 
or a known function of any of the other variables , equations (l) , (3) > 
and (4) comprise a determinate set, and hence the velocity u can be 
found. Before proceeding to discussions of F and the solution of the 
equations, however, the term containing pu in equation (3) requires 
further discussion. The result of the approximate equation (l) is the 
relationship in equation (2), which determines the product pu if it is 
known at any x. One must keep in mind, however, that, in the equations 
for a compressible gas, the density is related to the temperature, pres- 
sure, and velocity through the energy equation and the equation of state, 
as well as the continuity equation and the momentum equations. (These 
relations are discussed in a later section.) If the gas flow is to 
become fully developed, the density will assume a certain profile governed 
by these equations. It is not to be expected that pu will have exactly 
the same distribution profile at x = 00^ where the flow becomes fully 
developed, as it happens to have at x = 0. It is therefore evident that 
pu is not strictly a function of r, but may also be a slowly varying 
function of x. These considerations, however, do not invalidate the 
use of equation (2) in finding the velocity and pressure from equations 
(3) and (^-) in the regions where du/dx is not zero for the following 
reasons. The only instance where pu occurs in the momentum equations 
is in the term pu(du/dx) . This term, therefore, is approximated by 
PoUo(du/dx) . At x = 0, the acceleration du/dx is large, but the term 
is exact because pu = p 0 Uo • As x increases from zero, pu remains 
nearly poUo as du/dx decreases. As x approaches infinity, the 
actual pu profile approaches PoqU^ but du/dx approaches zero and 
the term becomes very small. Thus it is seen that the term pu(du/dx) is 
important near x = 0 where pu * Po^O; but its exact value is not as 
important at larger values of x. 


Body Forces and Electromagnetic Phenomena 


In the application of the solution of equations (2), (3)> and (*0 
to the magnetogasdynamic problem, F is the axial electric body force 
induced by electric and magnetic fields applied perpendicular to the 
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flow direction. If the magnetic permeability of the fluid is assumed to 
be the same as that of empty space, the electric body force per unit 
volume is 


F = J X 1 (5) 

— ^ — > 

where H is the strength of the magnetic field and J, the current 

density, is given by Ohm's law: 

J = cr(E + q X h) (6) 

Figure 2 indicates schematically the electric and magnetic fields . A 
uniform magnetic field, H 0 , is applied in the z direction. An appropri- 
ate distribution of electrical potential is applied to the conducting 
walls of the pipe . The fields are considered to be applied starting at 
the pipe inlet, x = 0. The interaction of the electric and magnetic 
fields is governed by Maxwell ' s field equations . The fields are , in 
turn, coupled with the fluid velocity by Ohm's law. It will be assumed 
that the z component of H does not vary appreciably from the value 
H 0 and that the magnitudes of the components induced in the x and y 
directions are small. These assumptions require that the magnetic Reynolds 
number cru^a be low. This has been discussed by Rossow (ref. 10) and by 
Hains and Yoler (ref. 13) among many others. Equations ( 5 ) and ( 6 ) then 
give for the axial body force 

F = crH 0 (Ey - uHq) (7) 

— > 

The electric field E is due to the internal separation of charges 
(polarization) as well as to the effect of the charged boundaries . It 
may be helpful to look at equation ( 6 ) in the form 

q. X H = ^ J - E e - Ei ( 8 ) 

where Ee is the electric field due to external effects (charged bound- 
aries) and is due to separation of charges. Thus,, from equation ( 8 ), 

it is seen that the induced emf, q X H, has two effects: that of driving 

a current through the gas and that of separating charges (from discussion 
by Resler and Sears in ref. 4) . 

The problem at hand will now be considered from two points of view. 

In case I it will be hypothesized that the current density component Jy 
is constant (fig. 2(a)). (Note that the purpose of this hypothetical 
case is to give a simple form for the body force , and it will probably 
be difficult to achieve in the circular pipe flow. Resler and Sears in 
reference k have shown it to be nearly true for two-dimensional channel 
flow through crossed fields.) Then, for case I, the electric body force, 
from equation (7) j becomes 

Fj = aH 0 (Ey - uHq) = constant 
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where a 
value , or 


where 


and where 


is also assumed to he a scalar and to have a constant finite 



Bl 


J y H o _ ( % 

cru 0 H 0 2 VH 0 u 0 



= constant 






(9) 


V 


K 2 . gH o 2a2 


The symbol K is the Hartmann number (ref. l4) , whose square is the 
product of the Reynolds number., the magnetic Reynolds number, and the 
magnetic pressure number. In this case E y , the y component of the 
electric field, varies because u is a variable. Case II (fig* 2(b)) 
will be the situation where the configuration and design of the apparatus 
and the magnitudes of the fields are such that E y can be assumed to be 
nearly constant. This can conceivably be true for large applied electric 
fields and for appropriate conducting properties and construction of the 
pipe wall. Thus, for case II, it is hypothesized that E y is constant, 
and the electric body force per unit volume is a variable and has the form 


where 



B II “ 




> 




( 10 ) 


The Approximate Dimensionless Equations of Motion 


Equations (2), (3)> and (4) will now be studied for both case I and 
case II where the body force has the form given by equations ( 9 ) and (10) , 
respectively. The parameters may be made dimensionless as- follows: 


— u 

- ap 

p 0 u 0 a 


u = — 
u o 

p Too 

N 0 = 

0 P 

* 

XI 

11 

P = f 

H o 

— /— \ am pua M -- 

m(r) = -— = ^ = N 0 pu 

* 

r = £ 

F - 

yz _ aHo 2 a 2 


a 

P-u 0 
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Then equations (3) and (4), using equation (2), become 



( 11 ) 



( 12 ) 


where, for case I 


F = Fj = K 2 Bj 


( 13 a) 


and, for case II, 


F = Fjj = K 2 (B n - u) 


(13b) 


A consequence of the foregoing hypotheses and assumptions is that 
the energy and momentum equations are not coupled, so they can he treated 
separately, and the momentum equations are linear after substitution of 
the approximation pu(du/dx) ~ p 0 uo(du/chc) . 

SOLUTION OP THE EQUATIONS OF MOTION 
Case I (F = K 2 Bj = Constant) 

Fully developed flow .- As x approaches infinity the flow approaches 
the fully developed flow condition and Su/Sx approaches zero. If this 
condition is imposed on equations (ll) and (12), and if the dimensionless 
velocity for fully developed flow is denoted as u = u^(r ) , the resulting 
differential equation for case I is 


where ( ) indicates differentiation with respect to r, and where: 


jf ( ru oo’ ) = -^ A 


(1*0 



(15) 


P = PooM 


(16) 


-r=r = constant 
dx 


(17) 
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Equation (l6) results directly from equation (12), and equation (l7) 
is easily deduced by observing that all other parts of equation (l4) are 
functions only of r. The boundary conditions are: 

r = 0 , Uqq is finite (l8a) 

r = 1 , Uoo = 0 (l8b) 

The solution of equation (l4) , with the boundary conditions (l8a) and 
l8(b), is 

u^(r) = A(1 - r 2 ) . (19) 

where A is defined by equation (15) . 

Development of the flow from a uniform profile A new variable U 
may be defined by 

U(x,r) = Uoo(r) - u(x,r) (20) 

where the function uio(r) is the fully developed dimensionless velocity 
profile, and u(x,r) is the actual dimensionless velocity at a point in 
the pipe between the inlet, x = 0, and a high value of x where the flow 
is considered to be fully developed. Substitution of equation (20) into 
equations (ll) and (12), using (l3a) , gives 


m 


du 

dx 


-Ft 


dp 4 d 2 U 1 dU d 2 U /I 


+ — + — 

dx 3 bx= 


+ ~ 0 
r dr dr 2 


- u. 


dp 


d^ 


1 

dr 3 dx dr 


= 0 



( 21 ) 

( 22 ) 


If equation (2l) is differentiated with respect to r, equation (22) is 
differentiated with respect to x, and the resulting equations are com- 
bined, the result is 


afu + i si - 1 M + JfSL - 1 (s *2) = o 

dr 3 r dr 2 r 2 dr dx 2 dr dr \ dx' 

where use has been made of the fact that, from equation (l4), 
(uco' 1 + (l/rOux,' ) is a constant . 


(23) 


In the present problem, the velocity at the pipe inlet is uniform, 
that is, u(0,r) =1 (0 < r < l) . The density profile at the inlet is 
also considered to be uniform, so that p"(0,r) = 1 and m(r) ~ mo(r) = N 0 , 
a constant, (0 < r < l) . The boundary conditions on u are thus: 
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A 


x = 0 , u = 1 , 


(0 < r < 1) 


x = “ ; u = UooCr) 

r = 0 y u is finite 




r = 1 , u = 0 


, J 

The corresponding conditions on U are, therefore: 


x = 0 , 

U = Uoo( r ) - 1 > 

(0 < r < 1) 

x = 00 , 

0 

II 


r = 0 , 

U is finite 


•s 

H 

II 

\U 

O 

II 



( 24 ) 


(25a) 

( 25 b) 

(25c) 

( 25 d) 


Particular solutions to equation ( 23 ) may he assumed to have the 
product form, X(x)R(r) . The complete solution is the linear combination 
of all possible solutions. Substitution of the product XR into equation 
(23) then gives the result 


U(x,r) = y A n e 6llX R n (r) 


n=i 


(26) 


where 


-n = § N 0 


1 - / 1 + 


2A 


•n 


W r 


( 27 ) 


and where A n 2 is the separation constant. (The positive sign on the 
radical has been eliminated by the boundary condition (25b).) The func- 
tions R n are the solutions of 


+ (j R* j + *n 2 R’ = 0 


(28) 


or 


R*' + = R' + A n 2 R = 


a 


(29) 


where a is a constant of integration. The eigenvalues, A n , are to be 
determined. The boundary conditions on R are obtained from equation (26) 
with conditions (25c) and ( 25 d) . Thus 
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and 


U(x,l) 


00 

:^A n e 6nX R n (l) = 0 
n=i 


U(x,0) 


= ) A n e €nX R n (0) 


n=i 


is finite 


Since "both these statements must be true for all values of x, it must be 

true that A 


R n (l) = 0 (30a) 

and that 

R n (0) be finite (30b) 

The solution to equation (29) is 

Rn = c iJo(^n r ) + C2Yo(^n r ) + ( 31 ) 

An 

where J Q and Y 0 are Bessel functions of the first and second kinds of 
order zero, and Ci and c 2 are arbitrary constants. 

Equation (31) involves three constants of integration, but only two 
direct conditions on R n (r) and its derivatives are available for the 
evaluation of the constants. These conditions, equations (30a) and (30b), 
will be used to evaluate ci and c 2 . The remaining constant, a, will then 
be chosen so that other conditions of the problem are satisfied. This can 
be done indirectly by considering all possible values for a and deter- 
mining the applicability of the results for each choice . The two classes 
of possible values for a that will be considered are: (l) a = 0 and 

(2) a f 0. 

Consider first the possible case a = 0. Application of conditions 
(30a) and (30b), along with the condition a = 0, to equation (31) gives 
the results: 


c 2 = 0 


and 


J 0 (An) =v0 


(32) 

(33) 


Equation (2 6 ) may then be written 


ON VO UO 


< OOCT\\0 


15 


( 34 ) 


U(x,r) = 2 j A n e £nX J 0 (A n r) 

n=i 

where the eigenvalues, A n , are the roots of equation (33) • 

Application of boundary condition (25a) to equation (34) leads to 
the result 

00 

(35) 


£( r ) ~ ^AnJp(A n r) 


n=i 


where 


6(r) = ■Uoo(r) - 1 , (0 < r < 1) (36) 

The functions Jo(^r?) are orthogonal with respect to r (cf. ref. 15) . 
Therefore equation (35) represents the expansion of the function |(r) in 
a series of orthogonal eigenfunctions (see again ref. 15) * and A n may 
be determined by the conventional method. Thus 

f 1 £(r)J 0 (A n r)r dr 

An = (37) 

f Q L J 0 ( 4 n r ) ] r dr 

Using equations (19) * (36), and (37)* one obtains 

. 2(kA/^ - 1 ) 


Thus the constant ci has been determined indirectly. It has been absorbed 
in A n * which is determined by equations (33) and (38) • 


The dimensionless velocity u(x,r) is now completely determined, for 
the constant a equal to zero, from equations (l 9), (20), and (34), along 
with equations (15)* (27)* (33)* and (38). 

One must consider now the possibility that a has a value other than 
zero. It will simply be noted here that it is shown in the appendix that 
the solutions obtained by letting a have any value other than zero are 
not possible . It is therefore established that 
\ 

a = 0 (39) 


The pressure may now be found as follows : Integration of equation 

(22) gives 


(40) 


P = 


5 1 + <*» 


1 6 


where G(x) is to be determined. Then, if equation (4o) is differentiated 
with respect to x and substituted into equation (2l), and use is made 
of equations (l4), (l5)> (26), (27 ) , and (29), the following result is 
obtained: 

dii 

G'(x)=-^ (Li) 


Integration of equation (4l) and its substitution, along with equation 
(34) , into equation (40) gives an expression which determines the pressure 
at any point to within a constant of integration. It would seem desirable, 
then, to find the constant of integration using the condition that, at 
x = 0 and r = 0, p = p Q . The following expression would result: 


P - Po 



— A e 
3 ■ H n fe n 


e n x n 


|(NjT) 



(te) 


Consider the case where dp^/dx = 0 (this will later be shown to be the 
only possible case in fully developed flow). Equation (42), with the 
condition x -► <*>, p -» p , then gives 



i—j j 
n=i 


(43) 


However, it is found that the above two equations are not convergent. But 
if the local press -are is related to the pressure at x = oo, for example , 
by adding equations (42) and (43), a result is obtained which converges 
for all x except x = 0, because of the presence of the exponential 
term: 


P - Poo = ' 5 A n € n e€nXj o( A n r ) W 

n=i 

The pressure difference, (p Q - p^) , can then be found by calculating 
(p “ Poo) ^ rom equation (44) for small values of x and observing the 
limit approached as x approaches zero. 
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Case II [F = K 2 (Bu - u) ] 


Fully developed flow.- For this case, equations (ll) , (12) , and (l3h) 
result in 


where 



C - Bjj 


]_ ^00 

K 2 dx 


( 45 ) 

( 46 ) 


Poo = pJ x ) 



constant 


( 47 ) 


With the boundary conditions (l8a) and (l8b) , the solution to the fully 
developed flow equation, (45), is 


Uco(r) = C 


1 


Ip(Kf) ~ 

Io(K)_ 


( 48 ) 


where I 0 is the modified Bessel function of the first hind of order zero 


Development of the flow from a unif orm profile . - As in case I , a new 

variable U(x,r) is defined by equation (20) . For this case u^(r) is 
given by equation ( 48 ) . Then, by an exactly analogous procedure as before 
equations (ll) , (12) , and (l3b) give the partial differential equation: 


^ + 1 ^ 2 U 1 6u 


. ^ — . K 2 d^U 
Sr 3 r 8r 2 1,2 8r 8r 8x 2 8r 


or \ ox- 7 


As in case I, the velocity and density profiles at the pipe inlet are 
uniform so that 


u(0,r) = 1 , (0 < r < 1) 


and 


m(r) = N 0 , a constant, (0 < r < l) 


The boundary conditions are then identical to equations (24) and (25) . 
The solution to equation (49) is found by the same procedure as used in 
case I to be 


l 8 


U = ^A n e enX R n (F) 
n=i 


where e n 


is now given by 



(50) 


(51) 


and the functions R n are the solutions of an equation identical to 
equation ( 29 ) • The equations in the subsequent development are the same 
as equations (30) through ( 37 )* Equation ( 48 ) is then used along with 
equations (36) and (37) to obtain 


A n = 


2 

AnJi(^n) 


C K_ 

A n 2 + K £ 



( 52 ) 


In the development of equation (52), the integral jj^ 1 1 0 (Kr) J 0 (A n r)r dr, 

which is similar to the Lommel integrals, was encountered. Its integration 
was handled in a manner following that used by McLachlan in reference 1 6 

( pp . 9^-95) • 

The dimensionless velocity u(x,r) is thus determined for case II by 
equations (20), ( 3 * 0 * and ( 48 ), along with equations (33), ( 46 ), (51), and 
(52) . The impossibility of a having any value other than zero is shown 
in the appendix. 

Following the same procedure used for case I, the pressure is found 
to be determined by an equation identical to equation ( 44 ), where A n 
and e n , of course, correspond to case II. 


ENERGY CONSIDERATIONS 

Equations of State and Energy in Fully Developed Flow 


Fully developed flow is considered to be a parallel flow in which 
the velocity profile does not change with streamwise distance . Maslen, 
in his study of fully developed gas flow in reference 12, states that, 
for a gas, none of the properties must be allowed to vary with axial 
distance. He obtains the result, then, that the velocity is only in the 
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x direction; that is, the flow is parallel. However , in the present 
study fully developed flow is defined exclusively in terms of velocity: 
q = iu where Su/Sx = 0, and where i is the unit vector in the x direc- 
tion. The author prefers to use this "nonaccelerating parallel flow” 
definition of fully developed flow, rather than a statement that nothing 
depends on x, in order to keep the defining statement identical to that 
used for incompressible fully developed flow, in which flow properties 
can vary with x. (Note that incompressible fully developed flow need 
not be defined to be parallel because it will be parallel if du/dx = 0 
to satisfy conservation of mass .) With this definition the result will be 
obtained that the flow properties cannot vary with x in fully developed 
laminar compressible viscous pipe flow. 


In magneto-fluid dynamics the Joule heating term must be included in 
the energy balance on a fluid particle . The energy equation to be used 
is derived from a classical macroscopic point of view. The fluid is 
treated as a single -component system, and therefore no consideration of 
electrons, ions, or neutral particles is introduced (ref. IT)- If the 
definition of fully developed flow (q = iu and du/dx = 0), which leads to 
u = u 00 (r) and v r = 0 is used along with the results of the momentum and 
continuity considerations, p = p 00 (x) (i.e., p is at most a function only 
of x) and dp^/dx = constant, then the energy equation in fully developed 
laminar compressible viscous pipe flow of a conducting fluid (assuming k, 
p, and Cp to be constant) is 


ST dp °o ^ V 

pucp - — = u + k 


chc 


dx 


d 2 T 1 dT d 2 T 

' 2 r dr dx 2 . 


+ H 


du 

,dr. 


+ ft 


where u = u^r) and where ft is the Joule heating term: 


V 


The equation of state is 


p = pRT = 7 ■ ■ ■ 1 pcpT 


(53) 


( 5 ^) 

(55) 


Because large temperature differences can exist in the flow in some cases, 
the fluid properties may vary considerably. The properties k, |i, Cp, 
and c r have been assumed to be constant in order to conveniently study the 
flow. The gas constant R is assumed constant, which implies that the 
number of particles per unit mass does not vary. 0hm T s law, relating the 
electric field to the magnetic field and velocity (assuming H 2 = H 0 )> is 

Jy = cr(Ey - uH 0 ) (56) 
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The parameters may he made dimensionless as before, with the 
additional parameters: 


- c p T 

T = ~2 3 

V-O 


M-Cp 


_ a. 
ft = — 




ft 


Equations (53) and. (55) then reduce to: 

TT _ dp^ 1 1 df S 2 t 

^0^00*^00 — ^oo l -n 2 — -x x 2 

ox dx Pr \ox r or or 


iu *f + 5 

dr/ 


(57) 


and 


where 


Poo (x) = N q ^ 


7 - - - 

PooT 


° ^ ^ 


(58) 


( 59 ) 


Two cases have been considered in previous sections where it has been 
hypothesized that only axial body forces exist. The current density J 
is assumed to be only Jy, and cr is assumed constant. In the first 
case, Jy is considered to be a constant. Thus, for case I, 


0 = ft j = constant 


(60) 


In case II, Ey is hypothesized to be constant. In that case, then 


ft = ftjj = K 2 (Bji - Uoo ) 2 


(6l) 


where B TT = Ey/H 0 u 0 . Also considered here will be case III, conventional 

gas pipe flow, that is, the case where a = 0 or K 2 = 0, with the result 
here that ft = 0 . 


Temperature and Density in Fully Developed Flow 
Corresponding to Case I 

For case I, discussed above, combination of equations (57) > ( 58 ), and 
( 60 ) and use of the fact that (p^u^o) is a function only of r give 
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3 

} 


i <3j?oo y Poo( x ) A , 


Uco — = = 


7-1 dx 7 - 1 Pr I 0 
where % is a constant and 


= Voo + Voo 


+ (Uoo') 2 + % 


“ (r) * pjr) 


(62) 

( 63 ) 


There are only two possible ways in which equation (62) can be satisfied. 
It is known from the momentum equations that dp^/dx: is a constant . 

Therefore, the only function of x in the equations is p^x) • IP is 
therefore concluded that either p^, is a constant or that 


1 » 


Pr N r 


- Voo + Voo 


= 0 


If the latter choice is considered, equation (62) becomes 


(uoo* ) 2 


7 - 1 J dx 


1 A dpoo 

T=r Uoo + ^1 = 0 


It can be seen immediately that this equation is not compatible with the 
momentum equation. Its integration, in fact, yields an impossible result. 
Therefore this case is not possible, and it must then be true that p^ 
is a constant. The equation of state, (58) , may then be rewritten, with 
p^ = constant. 


~ 7 - 1 NoPjF) 

to show that T can be a function only of r, not of x, that is, 

T = T M (r) (65) 

Therefore T^ is required to be constant . 

The density and temperature profiles in fully developed flow may now 
be found from equation (62) as follows: 

v «' ' + | ▼«,' = %(?) (66) 


( 61 .) 


where 
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= (^T 2 ) ^ Cfto’ ) 2 + %] (67) 

' ' ' ^00 

Using the result for for case I from the momentum equation: 

Uoo = ^ K 2 Bj(l - r 2 ) (68) 

and the fact that 

Si ’ k2 (i^' s ) ~ kV m 

equation (67) becomes 

* * Hr) (rsf 1 ) ^ d + gt) (TO) ' 

Equation (66) may then be integrated, using the conditions 


Voo(o) finite 


Vco(l) = N 0 ( 2 -^) ^ 


7 


Poo 


(71) 

(72) 


(where eq. (72) is found from the equation of state, ( 64 )) to give the 
fully developed flow density distribution: 


_1_ 

Poo 


7 - 1 No 
7 Poo 


T, 


«V + 4 * K B i 


4t , 2 fl - r 4 1 


f2 


16 


K* 


( 73 ) 


Equation ( 7 l) may then be combined with equation ( 64 ) to give a relation 
for the temperature profile in the fully developed flow: 


Too _ _ 1 


Pr 


K^Bi 


2 fl - r* 1 - r2 


16 


K“ 


( 74 ) 


Another variable of interest, the dimensionless mass-flow distribu- 
tion, which was discussed earlier and which will be of use in a later 
section, is 


cn o\vo 
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moo(r) = NqPooUoo = 


7 


- 1 Pr K 2 Bj. 


1 - ? 2 


4t 0 




i (l - r 4 ) + i (1 - 


Pr(K 2 Bj) 2 1 6 


K" 


r 2 ) 


(75) 


Temperature and Density in Fully Developed Flow 
Corresponding to Case II 


For case II, equations (57), ( 58 )? and- (6l) result in 

(= Voo' + Voo M ^) + (uoo*) 2 + K 2 (B n - u^) 2 


1 - ^00 = 7 Pqq(x) 

-l^dx 7 - 1 Pr N 0 


( 76 ) 


From the momentum equations, dp^/dx is a constant. Therefore the only 
possible function of x in equation ( 76 ) is p^x) . It again is con- 
cluded that either p^x) is a constant or (l/Pr N 0 ) [ (l/r)v 00 t + Voo' ' ] = 0. 
The latter choice implies no radial heat transfer; that is, either k = 0 
or T f T(r) . With this consideration, equation ( 76 ) becomes 


foo ’) 2 + k 2 (bh - ^) 2 


1 Vdx 


% K - 0 


( 77 ) 


But the solution for the fully developed flow velocity found from the 
momentum equation is 


Uoo = C 


1 - 


Ip(Kx) 

Iq(K). 


( 48 ) 


where 


n _ *d _ _L dp °° 
c " Bn p ax 


(46) 


For finite Hartmann number, equation (77) is contradictory to this solution 
(eq. (48)). However, let us consider for a moment the condition that 
Hartmann number approach infinity. Noting that 


d-Poc _ K dp OT 

dx ou 0 H 0 2 dx 


(78) 
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and that u^ and Ua,' approach C and zero, respectively, in (0 < r < l) 
as K approaches infinity, it is seen that equations (77) and (48) are 
then compatible in the region (0 < r < l) if 


dp oEyH 0 
dx / 


(79) 


so that 


7 - 1 „ 

u<x> - — “ — Bjj , 


f (0 < r < l) 

Ik -» oo 


(8o) 


Because p cannot change in the x direction, equation (79) implies also 
that the temperature have a constant positive gradient in the x direc- 
tion. This solution is not valid at the wall, that is, in the boundary 
layer, but it may be used in comparing the results of the present study 
with inviscid solutions in the case of no heat transfer. The flow 
described by equations (79) and (8o) cannot, of course, be considered to 
be a fully developed viscous flow because it requires that p = 0 . 


For case II, as for case I, therefore, the condition for fully 
developed compressible flow with nonzero viscosity is that p TO _be a 
constant and that dp^/dx =0. It is then also required that T = Too(r) 
and thus that be a constant . For the purpose of finding the density 

and temperature profiles, equation ( 76 ) may be written 


3" + - Voo' = T)jj(r) 


(81) 


where 

411 = C — ~~ -°' [(uoo' ) 2 + K 2 (Bh - vO 2 ] (82) 

' ' ' ^00 

and where, from equation (48), 


Uoo = B II 


1 


Ip(Kr) " 

I 0 (K)- 


( 83 ) 


After substitution of equation (83) and the relation 


Uoo' - "BjjK 


Ii(Kr) 

Io(K) 
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equation ( 82 ) becomes 

TluC?) = 


A - 7^ Pr N 0 

PlI K ] 

W ) p„ 

Llo(K)j 


[l! 2 (Kr) + I 0 2 (Kr) ] 


(84) 


Equation ( 8 l) may then be integrated, using conditions identical to 
equations ( 7 l) and ( 72 ), to give 


1_ 

P™ 


7 


1 “o % + ^If 

'- L °°w + 2 


1 - 


I 0 g (Kr) 

i 0 1 2 (k). 


(85) 


Use of equation (64) then gives the temperature profile for fully developed 
flow, case II: 


Too " 2 Pr Bjj 2 

The distribution of m^r) = NqP^u^, is 


7 PooBu 


moo(r) = 


7-1 


1 - 


_ Ip 2 (Kr) 

Io 2 (k) J 

Io(Kf) 
Io(K) 


(86) 


^ + 2 Pr Bl1 " 


1 - 


I 0 2 (Kf) 

Io 2 (K) J 


(87) 


Case III (a = 0) 


It is of interest to investigate the limiting case where the electri- 
cal conductivity approaches zero, that is, the case of conventional 
gasdynamic pipe flow. If it is supposed that such a flow can become 
fully developed, the result from the momentum equations is 


- K- W (1 - ?2) 


d Pc 


( 88 ) 


where dp /dx is a constant. The energy equation and equation of state, 
( 57 ) and ( 58 ) , for this case (where ft = 0) reduce to 

1 - d Poo 7 Poo(x) /l , . ft ^ , tN a ( q ^ 

Uoo ~ 7 = = r - w ■ (=■ Vo,' + Voo" ) + (uco* ) (89) 


7-1 


dx 7 - 1 Pr N 0 


2 6 


where v^, = l/p^. Because p^x) is the only possible function of x 
in equation (89) , it must be true that either p^, is a constant or 
(l/Pr N 0 ) [(l/r)^' + Vqo' '] = 0 . Substitution of the latter consideration 
into equation (89) results in 



which is contradictory to the momentum equation, besides the fact that 
it implies that Uoo' is imaginary if the necessary negative pressure 
gradient is assumed to exist . The only other possibility for fully 
developed flow, then, is p,*, = constant. But then dpoo/dx = 0 and 
equation (88) then states that there is no flow. The conclusion is that 
it is not possible to have fully developed laminar flow in a pipe if the 
fluid is compressible unless the fluid particles are acted upon by body 
forces, that is, forces other than the surface forces due to pressure 
and friction. This is perhaps a surprising result because it is known 
that air flowing at low speed in a pipe does, at least very nearly, 
approach fully developed flow. But at low speeds viscous dissipation of 
energy is small and the flow is nearly incompressible. The momentum 
equations can then be satisfied without recourse to the energy equation; 
the pressure gradient is not required to be zero, (in the case where the 
fluid is not exactly incompressible, but small density changes with 
temperature, i.e., thermal expansion and contraction, are admitted, Maslen 
(ref. 12) has obtained solutions using an equation of state which is 
independent of pressure and using momentum and energy equations which 
involve a constant pressure gradient.) At high speeds, in addition to 
the fact that the density must be allowed to vary, there is greater 
viscous dissipation of energy, which fact lends to acceptance on physical 
grounds of the result found above, that a fully developed flow is not 
attained unless it is supported by a body force. 


MASS -FLOW COMPARISON 


In the section "Equations of Motion" physical arguments were given 
in support of the use of the hypothesis that v r = 0 and the consequent 
use of the function p Q u 0 (du/dx) as an approximation for the term pu(8u/dx) 
in the momentum equation. The approximation is, of course, exact at x = 0 
because (pu) x _ 0 = p 0 u Q . It is also true as x approaches infinity because 

8u/8x approaches zero. Use of the approximation at values of x near 
zero is justified because, although u changes rapidly from Uq, the 
product pu remains nearly at p 0 u Q . At large values of x the approxi- 
mation is good because 8u/8x becomes smaller and hence the term pu(8u/8x) 
in the differential equation becomes less important . The approximation is 
least accurate at intermediate axial locations, that is, at values of x 


CT\VO uo > 
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where the center -line velocity is approximately intermediate between 
u Q and and the product pu is approximately intermediate between 
p 0 u 0 and Poo^oo* One can obtain, therefore , an indication of the accuracy 
by comparing the mass flow at x = °o with the mass flow at x = 0, that 
is, by calculating the ratio m^/m 0 versus r. One can then consider 
the value of S/Sq intermediate between S»/mo and 1.0 to be a good 
indication of the accuracy of the approximation at the x location where 
it is least accurate. For this calculation equations (75) and ( 87 ) are 
used to find for cases I and II, respectively, and m^ is equal 

to N 0 . 


First it must be required that over- all conservation of mass be 
satisfied. The respective expressions for m^ contain the parameters 
Poq and T^ which must have definite values depending upon the state of 

the flow at the pipe entry. Thus the various parameters must be related 
so that they satisfy the equation 


or 



( 90 ) 


Case I (J v = Constant) 


For case I, the result of substituting equation (75) into (90) 
requires that the parameters be related by 

1+S 3 

(l + &3 + 6 4 )( 5 3 - 64)1 5 


6-1 = In 


bp + bg + TT 


(l + S3 - S 4 )(S 3 + s 4 )_ 


( 91 ) 


where 


bp - 


81 = 

32^ 


EU ( Pr )(K 2 Bj) 


8 P„ 


32 ^oo v /To 


Pt(K 2 Bj) 2 


7 — Pr )7M 0 2 


(K 2 B!) 2 


(92) 

(93) 
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(94) 




63 + 263 + 2£>2 + 1 


(95) 


As will be discussed later, the results of calculations demonstrate the 
fact that 



K 2 Bj 
P °° = 3N 0 


Combination of equations (92) and (96) yields the relation 


(96) 


K 2 Bj = 


3 Np 2 861 

7Mq 2 8&! + 3No 2 Pr(7 - l)/7 


(97) 


From equations (9l) and (97) a definite relationship between K 2 and K 2 Bj 
can be seen for given values of N 0 , [(7 - l)/7]Pr, 7Mo 2 ; and T^/To- 
One may note that, for those cases where 

K 2 Bj « 3W 0 2 /7Mo 2 

equation (97) may be approximated by 

o 8& n 

(CT 1 Pr )( 7M ° 2 ) 

and the relationship between K 2 Bj and K 2 would then depend only on 
[(7 - l)/7]Pr(7Mo 2 ) = (d/k) ( u 0 2 /T 0 ) and on To^/To* This is true, for 
example, if N 0 > 100, K 2 Bp < 1000 and 7 and Mq are of order unity. 

The calculation can be made as follows: For various given values 

of K 2 , No, [(7 - l)/7]Pr, and 7Mo 2 , a number of arbitrary values of 
&2 are specified. Corresponding to each of these values, then, is 
known from equation (91) • The values of K 2 Bj are then found from 
equation (97) , and the appropriate values of T^/To are obtained from 
equation (93) • 

For chosen sets of the input parameters which satisfy equation ( 90 ), 
as calculated above, one may then calculate the quantity ma/mo as a 
function of r. Using equation (75) > one obtains 
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mo 



1 - r 2 

6 a + | (1 - r 4 ) + 6 3 (1 - r 2 ) 


(98) 


Case II (Ey = Constant) 

For case II, if equation (87) is substituted into equation (90) , the 
required relationship among the parameters is found to be 


where 



(l - £)r dr 

b 6 + 1 “ 


and where 


, _ Iq(KF) 
5 Io(K) 


85 = 


N 0 



Pr Bqj 



&6 


2Too w 
Pr Bjj 2 




(99) 

(100) 

( 101 ) 
( 102 ) 


In this case the results of calculations } to be discussed later , reveal 
the fact that 


_ k2 ( e ii - h 

7Mq 2 " P °° 3N 0 


(103) 


Equations (lOl) 


and (103) may be combined to give an expression for 


Bji: 


B II “ 



4S5K 2 + 3N 0 2 Pr 



(104) 
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For this case the calculation can he accomplished in a manner 
similar to that used in case I. For various given values of K 2 , N 0 , 
[Pr(/ - l)/7], and /Mo 2 , a number of arbitrary values of &e are 

specified. Corresponding values of 5s are found from equation (99) 

(by numerical, integration, for example), and then Bqj is obtained from 
equation (lO^) . The appropriate values of To^/Tq are obtained from 
equation (102) . 

Corresponding to the sets of input parameters which satisfy equation 
(90), as calculated above, moo/5^ may then be calculated. For this 
purpose equation ( 87 ) may be written: 

I • (4) G. t : - p) < iq 5> 

where £ is given "by equation (lOO) . 


NUMERICAL EXAMPLES 


The equations containing the infinite series solutions for the 
dimensionless velocity u and the dimensionless pressure (p - p^) at 
all x and 7 were programmed for evaluation on an IBM type 70k electronic 
data processing machine . Also calculated were the fully developed velocity 
and temperature profiles . The equations for both case I and case II were 
programmed . 

A subroutine employing Newton* s method as the iteration procedure 
was used to evaluate the eigenvalues, A n , as the successive roots of 
equation (33) • An approximate value was found graphically for the first 
eigenvalue, Ajl. This approximate value was then used as a trial, or 
starting, value in the iteration procedure to solve for the actual value 
of Ajl. The trial value for the iteration in the vicinity of each 
succeeding eigenvalue, A n , was found by adding the number it to the 
preceding eigenvalue, A . 

The calculations were made for values of the Reynolds number N 0 
of 1000, 750, 500; and 250, values of K 2 of 1000, 500; 250, 100, 50; 
and 10, values of Bj = (Ey/H 0 u 0 ) - (u/u 0 ) = JyH 0 /ou 0 H 0 2 of 1.0 and 0.1 
and values of Bjj = Ey/H 0 u 0 of 100, 10, and 2. An appropriate range 
of x for each set of input parameters was used so that the asymptotic 
approach of the velocity distribution to the fully developed flow profile 
could be observed and so that the pressures near x = 0 could be 
determined . 


00 on vo 


o\vo U) 


31 


DISCUSSION OF RESULTS 
Fully Developed Flow 

The fully developed flow velocity, toward which the flow is assumed 
to accelerate starting from the inlet velocity, depends on K 2 and Bj 
for case I and on K 2 and Bj j for case II . Plotted in figure 3 are the 
velocity profiles referenced to the value of the velocity on the pipe 
axis. As presented in this figure, the profile is independent of all 
parameters for case I. In this case the electromagnetic body force, which 
is balanced by the viscous forces in the fully developed flow, is constant 
over the cross section. The result, a parabolic profile, is the same as 
A for the Hagen-Poiseuille incompressible pipe flow in which a constant 

pressure gradient maintains the flow. The velocity profiles depend on 
the Hartmann number K in case II. In this case the body force varies 
over the cross section. It is greater near the wall and thus tends to 
flatten the profile. As K approaches infinity, the body force approaches 
zero everywhere in the fully developed flow except in a boundary layer 
near the wall. Shown in figure k are curves of the fully developed flow 
center-line velocity versus K 2 for various values of the parameters 
Bj and Bjj. In case I, where J y is constant and the body force is 
constant and thus has a nonzero value in fully developed flow, the center - 
line velocity is directly proportional to K 2 as well as Bj. (The 
product K 2 Bj is actually proportional to the ratio of the constant 
electric body force to the inlet viscous force on a fluid particle; that 
is, K 2 Bj = F/(puo/a 2 )*) Thus, for this case, as the magnetic field is 
increased, the maximum velocity is also increased because the body force, 
which maintains a constant value throughout the entire flow, is increased. 
In case II the flow -approaches the fully developed condition as the body 
force approaches a value proportional to K 2 Bji/Io(k) on the center line. 
Thus, for values of K near zero, the body force on the center line 
becomes proportional to K 2 Bjj. But as K approaches infinity, the body 
force approaches zero. The fully developed flow center-line velocity is 
given by Uoo/uo = Bn[l - l/l 0 (K)]. Thus, as K nears zero, the fully 
developed flow velocity approaches zero. But as K approaches infinity, 
Uco/uo approaches Bjp (except very near the wall) and therefore becomes 
essentially independent of K. 

The results for case II may be compared with those of Resler and 
Sears in their study of nonviscous one -dimensional magnetogasdynamic flow 
in a constant-area channel with no heat transfer (ref. 5). Resler and 
Sears distinguished essentially three basic types of flow (with respect 
to velocity): (l) choked flow, (2) smooth continuous acceleration (or 

deceleration), and (3) approach to an asymptotic state. When the viscosity 
is not zero, the flow properties (velocity, temperature, etc.) vary over 
the cross section, and the conditions which lead to choking as the velocity 
changes with x would therefore not be present over the entire cross 
section at any axial location. Hence the effects of choking are not 
evident in the present study and choked flow is not considered. Neither 
have the cases for which Resler and Sears found the flow«to accelerate 
indefinitely with distance along the channel been investigated here . 
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Only those cases for which the flow velocity approaches the fully 
developed, or nonaccelerating unidirectional, flow state have been 
considered. The boundary condition was imposed that the flow asymptoti- 
cally approach this state. For inviscid flow with no heat transfer, 

Resler and Sears found the asymptotic flow velocity to have the value 
E/H in some cases, while in others it was limited to [(/ - l)/7]E/H. 

If the effects of viscosity are small enough, the velocity and temperature 
profiles will be uniform, and hence one -dimensional, except in the bound- 
ary layer at the wall. All effects of viscosity and thermal conductivity 
will be confined to the boundary layer, and the flow outside the boundary 
layer will be governed by the inviscid flow equations of Resler and Sears . 
Thus, if u 0 < [(7 - 1)/7]E/H < E/H, as the velocity outside the boundary 
layer approaches the value [(7 - l)/7]E/H, the flow will indeed cease to 
accelerate and hence the velocity will become asymptotic to [(7 - l)/7]E/H. 
In the present analysis (case II), if the viscosity is assumed to approach 
zero,, then the Hartmann number becomes infinite and the flow does become 
essentially one -dimensional, except for the boundary layer near the wall 
(see fig. k) , and both of the asymptotic flow cases indicated by Resler 
and Sears (u -> E/H and u [(7 ■ l)/7]E/H) are shown to be possible 
solutions of the fully developed flow equations, under the proper respec- 
tive conditions, in the limit as K approaches infinity (see eqs . ( 80 ) 
and ( 83 )). However, if K is small enough that the velocity profile is 
not extremely flat, then the first-order differential equations of Resler 
and Sears do not apply, and the only possible asymptotic flow velocity 
is that given by equation ( 83 ) . The effects of heat conduction will not 
permit the flow to become fully developed as the velocity approaches 
1(7 - l)/7]Ey/H 0 , and hence the flow will continue to accelerate until 
the center-line velocity approaches asymptotically the value, 

(Ey/H 0 )[l " l/loOO] (provided the proper constant wall temperature is 
supplied). As noted, however, the limiting results for inviscid flow 
agree completely with the asymptotic flow cases discussed by Resler and 
Sears. 


The temperature distributions in fully developed flow are shown plot- 
ted in figure 5* The results, as presented here in dimensionless form and 
in terms of the various parameters as indicated, depend only on Hartmann 
number, K. The temperatures are referenced to the constant wall tempera- 
ture. It is seen that the temperature profiles are governed to a great 
extent by the nature of the distribution of Joule heating. The effects of 
the various parameters on the temperature profiles are very similar to 
their effects on the velocity profiles . In case I the Joule heating is 
everywhere constant. As K approaches zero, the temperature profiles 
approach parabolic shape. As K goes to infinity, the shape becomes pro- 
portional to [l- (r/a) 4 ]. Because the Joule heating is nowhere zero in the 
fully developed flow for this case, the temperature continues to increase 
as K increases. The Joule heating in case II varies over the cross 
section, being greater near the wall. As K approaches zero, the tempera- 
ture difference (T- T w ) approaches zero. As K approaches infinity, the 
Joule heating approaches zero everywhere except in the boundary layer, and 
thus the temperature profiles approach a uniform profile . 


CT\ VO UU 


0\V0U) > 


D 


33 


Velocity in the Accelerating Portion of the Flow 


The development of the velocity profile from the uniform profile 
at the pipe inlet is illustrated by several examples in figures 6 and 7 • 
Figure 6 gives examples corresponding to case I. For all cases here, the 
flow develops, starting with a thin boundary layer near the inlet, to 
the parabolic profile of fully developed flow (although the parabolic 
shape is not evident from the semilog plots of fig • 6) • The effects of 
the various parameters will best be observed by studying the center-line 
velocities, as the radial development is very similar for all cases. It 
will be noted, however, that figures 6(b) and 6(c) are identical because 
these results depend only on N 0 and on the product of K 2 and Bj in 
case I. Figure 7 illustrates case II. For high values of K, the bound- 
ary layer remains thin, but as K decreases, the boundary layer develops 
into a profile approaching a parabolic shape . 

The center-line velocity distribution along the axial coordinate in 
the accelerating portion of the flow is a function of the inlet Reynolds 
number N 0 and the product K 2 Bp for case I; in case II it depends on 
the values of N 0 , K 2 , and Bjj. Figure 8 shows the center-line velocity 
relative to the inlet velocity. It is noted here that for small values 
of K 2 Bj in case I, the velocity overshoots the fully developed flow 
velocity before approaching it asymptotically. In the accelerating flow, 
as in the fully developed flow, the body force in case I is constant; 
whereas in case II, it varies with axial distance as well as over the 
cross section. Thus the body force in case II is greater near the inlet 
than near the fully developed flow, where it approaches zero as K 
approaches infinity. Nevertheless, the center-line velocity change in 
the accelerating portion seems to behave nearly the same in case I as 
in case II; that is, the shapes of the curves are nearly the same. It 
is noted, however, that increasing K in case II causes faster acceler- 
ation toward fully developed flow, an effect which is not seen in case 
I because increasing K in case I also increases the fully developed 
flow velocity. In figure 9 the center-line velocity distributions are 
shown relative to the fully developed flow velocity. Figure 9(a) shows 
the effect of varying K 2 and Bj for case I with No having the value 
1000. As K 2 Bj approaches infinity the axial velocity distribution, in 
the form presented in figure 9(a), approaches a single curve. One will 
note that the curve for K 2 Bj = 1000 cannot be distinguished from that 
for K 2 Bj = 500* la figure 9(h) can be seen the effect of varying No 
in case I. Higher values of No cause the velocity to approach the 
asymptotic profile more slowly, that is, the fluid proceeds further along 
the pipe before the flow becomes fully developed. The effect of varying 
K 2 and Bji for case II is given in part (c). Increasing K 2 causes 
the flow to approach the fully developed condition more rapidly. Higher 
values of Bij cause the flow to reach the fully developed state more 
slowly because the fluid must accelerate to higher velocity. Part (d) 
shows that the effect of varying No in case II is essentially the same 
as for case I. 
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From an expanded plot similar to those used in figure 9 can be read 
the values of x*/a; that is, the values of x/a where the center-line 
flow velocity has reached a value within 1 percent of the fully developed 
flow velocity. It was found that the best type of plot for this purpose 
was one of log |1 - (u/‘u 00 ) r _ o | versus (l/N 0 )(x/a) for case I and 
log [1 - (u/uoo) r=0 ] versus (K 2 /N 0 )(x/a) for case II. On these plots the 
results were very nearly straight lines, except for portions of the curves 
for the "overshoot" cases of case I. When the results for x* in case I 
are plotted as indicated in part (a) of figure 10, they fall on a single 
curve for all values of No in the range 250 to 1000. It is seen then 
that x*/a is directly proportional to N 0 and, at high values of the 
product K 2 Bp, it becomes nearly independent of K 2 Bp. As K 2 Bi decreases 
to low values, near ^-.0, x*/a also decreases. The reason for this is 
that the fully developed flow velocity is also directly proportional to 
K 2 Bi, and is equal to (l/^)K 2 Bi on the pipe axis. Thus, as K 2 Bp 
decreases, the pipe length required to achieve the fully developed flow 
velocity on the axis decreases because the velocity change is decreased. 

In fact, when K 2 Bj is very small., the fluid decelerates to fully devel- 
oped flow, and, as K^p approaches zero, the value of x*/ a rises again 
and approaches infinity, because the flow must approach zero velocity 
because of the requirement of zero pressure gradient in the fully developed 
flow. For case II, plotting the results in the manner shown in figure 
10 (b) brings them to the same order of magnitude for the values of the 
input parameters used. It is seen that, for high Noj x*/a is nearly 
directly proportional to N 0 • It decreases with increase in K 2 because 
the body force is then higher near the inlet, but lower near the fully 
developed flow. In fact, for this case, x*/a is nearly inversely 
proportional to K 2 for fairly high values of K 2 . It increases with 
increase in Bpp (for Bpp greater than unity) because the final velocity 
approached is higher. As Bpp decreases from unity to zero, the value 
of x*/a must approach infinity because the fully developed flow velocity 
approaches zero and the fluid must then decelerate from the inlet velocity 
to zero. These results for both cases I and II support the conclusions 
reached in the study of the energy equation in a previous section; that 
is, compressible viscous flow cannot become fully developed unless it is 
supported by a body force. Thus as K 2 , Bp, and Bpp, the significant 
parameters involved in the strength of the body force, approach zero, 
the pipe length required to achieve fully developed flow goes to infinity. 

Values of x*/a in case I are, in general, higher than those in 
case II. In case I the flow takes longer to become fully developed 
because the body force is applied at a constant rate along the tube, and 
energy is being added to the flow through Joule heating at a constant 
rate. In case II the flow becomes fully developed more rapidly because 
both the body force and the energy added by Joule heating decrease as 
the velocity increases toward the fully developed flow value. 
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Pressure in the Accelerating Portion of .the Flow 


The center-line pressure distribution over the axial distance is 
shown in figure 11. The results are calculated from equation (44) using 
r = 0 and the fact that 

P - P„ = No 

When the results for case I are plotted in the manner illustrated in 
figure ll(a), they are found to he nearly independent of No for values 
of No in the range 250 to 1000. Results are shown for various values 
of the product K 2 Bj. It is seen that as K 2 Bj increases, the results 
as plotted here approach a single curve , and thus that the pressure 
distribution becomes proportional to K 2 Bj, or to the body force. For 
case II (fig. 11(b)) the results for No in the range 250 to 1000, K 2 
in the range 50 to 1000, and Bjj in the range 2 to 100 fall on a single 
curve • 

The pres sure* difference between the value at the pipe inlet and the 
fully developed flow value can be found by observing the value approached 
as x approaches zero in figure 11. It is noted that on both parts of 
figure 11, all curves approach the value 1/3 on the ordinate and therefore 
that the value of (po~Poo)/Po u o 2 is (i/3)K 2 Bi/No 2 for case I and is equal 
to (l/3)K 2 (Bjj - l)/No 2 in case II. Thus the difference between the 
pressure at the inlet and that in the fully developed flow can be calcu- 
lated directly. Although the relationships among the various parameters 
are obvious from these formulas, numerical results are presented in 
figure 12 to show the order of magnitude of the pressure difference for 
various values of the parameters. It is seen that the pressure differences 
required are in most cases very small compared to the inlet dynamic pres- 
sure. It must be noted that for a given set of input parameters. Ho, K, 
and Bp or Bjl, there is a definite pressure difference. Thus, if this 
pressure difference does not exist in the actual flow situation, then 
the indication is that the fully developed flow condition will not be 
attained . 



Results of the Mass-Flow Comparison 


For the purpose of obtaining an indication of the accuracy of using 
p 0 u 0 (du/dx) as an approximation for the term pu(du/c)x) in the momentum 
equation, the dimensionless mass-flow profile in fully developed flow, 
moo, may be compared with the dimensionless inlet mass flow, m 0 , as 
discussed in a previous section, "Mass-Flow Comparison." Shown in 
figures 13 (a) and 13 (b) , respectively, are examples of the required 
relationships among the parameters for case I and case II, which must 
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first be calculated (in order that over-all conservation of mass be 
satisfied) . Note that although the examples plotted in figure 13(a) were 
calculated for N 0 = 1000, the results are essentially independent of 
No for the ranges of parameters used if N 0 is greater than a number 
substantially less than 100. The results shown therefore also correspond 
to any combination of values for /Mo 2 and [Pr(/ - l)//] which have the 
same product as those listed. (See earlier text under "Mass-Flow 
Comparison. 11 ) 

Figures l4(a) and l4(b) show examples of versus r correspond- 

ing to several of the relationships given in figure 13 • Some specific 
examples to which the cases plotted in figure 1^ correspond are listed in 
tables I and II. For these typical example cases and from the considera- 
tions outlined under the section "Mass-Flow Comparison," one can determine 
for which sets of conditions the approximation under consideration is most 
valid. (Recall that the approximation is most correct for those cases 
where m/niQ is near unity.) Evidently the validity increases as one 
progresses through cases I through k for case I (see table i) and through 
cases 5 through 7 for case II (see table II ) . 


CONCLUDING REMARKS 


It has been found that for the occurrence of fully developed laminar 
flow of a compressible viscous fluid in a pipe, an axial body force must 
be applied. That is, the viscous forces must be balanced by fluid body 
forces rather than simply by forces due to a thermodynamic pressure 
difference. It is thus found that, as the effects producing the body 
forces are decreased to zero, the pipe length required for fully developed 
flow becomes infinite. 

It is also shown that, for the flow of a compressible viscous fluid 
to become fully developed, the wall temperature must approach a constant 
value, and it is seen that definite values of the pressure difference 
between the value at the inlet and the constant value in fully developed 
flow are required in order that the fully developed flow be smoothly 
obtained. The radial temperature distribution is also seen to approach 
an asymptotic profile. 

In the application of the solutions to the magnetogasdynamic pipe 
flow problem, in which the axial body force is the "electric body force," 
the flow results depend on the nature of the electric field and current 
density. Two special cases have been studied in this report: In case I 
the component of current density perpendicular to the velocity and the 
magnetic field is assumed constant. For this case the velocity profiles 
approach parabolic shape as the flow approaches the fully developed 
condition. The temperature in the fully developed flow is found to 
depend on the Hartmann number. As Hartmann number, K, approaches infinity, 
the temperatures become proportional to the fourth power of K, and as K 
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approaches zero, the radial temperature profile becomes parabolic and 
proportional to the square of K. The pipe length required to achieve 
a velocity within 1 percent of fully developed flow velocity is found to 
approach a value proportional to the inlet Reynolds number for high 
values of the Hartmann number. The pipe lengths are in the order of 
hundreds of radii for a Reynolds number of 1000. The difference between 
the pressure at the inlet and that in the fully developed flow is found 
to be very small in comparison to the inlet dynamic pressure, the ratio 
being proportional to the ratio of the constant electromagnetic body force 
to the inlet viscous, force on a fluid particle and inversely proportional 
to the square of the inlet Reynolds number . In case II the component of 
the electric field perpendicular to the velocity and the magnetic field 
is considered constant throughout the flow field. The velocity profiles 
for this case approach fully developed flow profiles which become flatter 
with increase in Hartmann number. As the viscosity approaches zero 
(Hartmann number approaches infinity) , the results approach those of the 
inviscid constant -area channel flow studied by Resler and Sears. In 
contrast to the two different values of the asymptotic velocity indicated 
by Resler and Sears for inviscid flow, it is shown that only one asymptotic 
flow velocity profile is possible for a given set of conditions when the 
effects of viscosity are significant, which profile corresponds to the 
upper limit indicated by Resler and Sears, modified by the Hartmann number. 
The fully developed flow temperature profiles for case II behave much like 
the velocity profiles. For high values of the Hartmann number, the entry 
pipe length is nearly proportional to the ratio of inlet Reynolds number 
to the square of the Hartmann number. It also increases as the ratio of 
the electric field to the product of the magnetic field and inlet velocity 
increases . The pipe lengths are in the order of 5 radii for Reynolds 
number of 1000 and Hartmann number equal to *Jl000 - In case II the 
difference between the pressure at the inlet and that in the fully 
developed flow is also found to be small, in general, compared to the 
inlet dynamic pressure • The ratio of the pressure difference to the 
inlet dynamic pressure is proportional to the ratio of electromagnetic 
body force at the pipe inlet to the inlet viscous force on a fluid 
particle and inversely proportional to the square of the inlet Reynolds 
number . 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Feb. 13, 1961 
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APPENDIX A 

EVALUATION OF THE CONSTANT a FOR THE FLOW 
DEVELOPMENT FROM A UNIFORM PROFILE 


The solutions to the equations of motion corresponding to the constant 
of integration, a, in equation (29) having the value zero have been given 
in the text of this report . The possibility of solutions for a having 
a value other than zero will be investigated in this appendix. 

If a is assumed to have a nonzero value, application of conditions 
(30a) and (30h) to equation (31) gives the result 


Rn 


-Q- Jp ( ^n r ) 

^n _ Jq('V) 


(Al) 


where the values of a and the eigenvalues A n are, as yet, undetermined. 
It will be seen that a need not be determined directly, since it auto- 
matically falls out of the equations. Equation (Al) may be substituted 
into equation (26) to give 


U(x,r) 


00 



n=i 


~’J 0 ( A n.r) 
_ Jo(^n) 


1 


(A2) 


where 


B n = 



(A3) 


Application of boundary condition (25a) to equation (A2) then leads to 
the result 


6(F) 


00 

) ^ bn J o ( ^n r ) 


(Ah) 


n=i 


where 


00 


!(r) = ujr) - 1 



n=i 


(0 < r < l) 


(A5) 


and where 


D n 


% 

Jo( A n) 


(a6) 
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If the functions J 0 (A n r) are orthogonal, then equation (A4) 
represents the expansion of the function £(r) in a series of orthogonal 
eigenfunctions. In order for Jo^n?) to be orthogonal, one of the three 
possible orthogonality conditions must be satisfied (see ref. 15 )* The 
most general of 'the three is 

pJ 0 (^n) + ^n^i(^n) = 0 {A’j) 

where 3 is a negative constant . It will be noted that the other two 
possible conditions, Jo(^n) = 0 an( i ^i(^n) = are merely special cases 
of equation (Aj) , that is, the respective cases where P = -oo and P = 0. 
The case where 3=0 must be ruled out because this would impose the 
condition, through equation (A2) , that dU(x,l)/Br = 0. This requires 
the restriction that at the wall du/dr = u^ 1 (l) at all x, which is 
unreasonable for the case of a uniform velocity profile at the inlet . 

The case 3 = -oo corresponds to a = 0, which case has been discussed 
in the text . 

The choice of various finite negative values of 3 in equation (A7) 
would supposedly result in different physical conditions of the problem. 
For example the difference in pressure at two points in the flow would be 
a function of the choice of 3 • Thus a complete set of eigenvalues are 
prescribed as the roots of equation (A7) with a given choice of 3; and 
the appropriate pressure distribution would result. 

The coefficients D n may then be determined by the conventional 
method: 

f 1 !(r) J 0 (A n F)r dr 

D n = ^ (A8) 

/ [Jo(A n r)] 2 r dr 

o 

Case I (F = K 2 Bj) 


For case I, using equations (l 9 ), (A 5 ), (A6 ) , and (A8) gives the 
result 


where 



f = f(A n ) = 


2Ji(A n ) 

Ant Jo(An) ] 2 + AntJ^An )] 2 


(A9) 


(A10) 


to 


and 


g = g(A n ) 


2Jp(A n ) 
AnJ i(^n) 


(All) 


If the sum over values of n from one to infinity is taken of equation 
(A9) , there results 



If r = 0 is substituted into equations (Ak) and (A5) , using equation (19), 
the following relation is found: 



00 


If the above two relations for 
results: 



are equated, the following expression 



(A12) 


Substitution of equation (A12) into equation (A9), making use of equations 
(A7) , (A10), and (All), gives the result 


_ 2A /g - 2 + 8(3 2 S 2 _ 4p\ 
n Aq 2 + p 2 \ 1 + 2PS X A n 2 / 


where 


and 


Sr 


V 1 fl/Jo(A n ) 
La LA n 2 + P 2 . 


00 



n=i 


l/Jo(An) 
-An 2 (A n 2 + P 2 )- 


% 

(A13) 


(Al4) 


(A15) 
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There are two indications that this solution which has "been found 
for a ^ 0 is not possible. First, it is found by computing Si and S 2 
that Si = -l/20 and that S 2 = (2 - 0)/80 2 . Therefore the first term 
inside the bracket in equation (A13) is indeterminate, that is, both the 
numerator and denominator of that term are identically zero for all 0 . 
The second indication is seen by comparing equation (A9) with equation 
(AI3) , using also equation (A6) . The entire right side of equation (A13) 
has the independent parameter A as a factor . Therefore , in order that 
equation (A13) be compatible with equation (A9), it must be true that 


either (l) ) B n = 1, or (2) 


B 


n 


must have a factor A. The first 


n=i |_n=i 

consideration cannot be true because, from equation (AI3) , B n has a 


factor A. Therefore, ) B n would also necessarily have a factor A. 


n=i 


This same point rules out the second consideration because if ) Bn has 

”00 ~] Z_J 

V -1 n=i 


a factor A, then 


B n ’ 1 


n=i 


cannot have a factor A. It is therefore 


shown for case I that the solutions corresponding to values of a unequal 
to zero are not possible . 


Case II [F = K 2 (Bjj - u) ] 

For case II, equations (48), (A 5 ) , (A6) , and (A8) lead to 



where f = f(A n ) is given by equation (A10) and 


(A16) 


h = h(A n ) = 1 



Ii(K) Jp(An) " 
Io(K) A n Ji(A n ) + _ 


(A17) 


If the sum over values of n 

00 


(Al6), a relation for 



from one to infinity is taken of equation 
results which may be equated to another 


n=i 
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relation for D n found by substituting r = 0 into equations (A4) and 


n=i 


(A 5 ) . A resulting expression is 


00 



n=i 


(A18) 


Substitution of equation (Al 8 ) into (Al 6 ) , using (A7 ) , (A 10 ), and (A17), 
results in 


( ,2f» 

~1 + 2ps 3 - l/lo(K) _ K 2 + (l/pjAn^iC^/loCK)" 

W 2 + p 2 / 

1 + 2PSi K 2 + A n 2 


(A19) 


where Si 


is given by equation (Al4) and where 


^ — t ( 1 /P) A n 2 KIi(K)/l 0 (K) 

V K 2 + A n 2 

/ Jo(^nH^n 2 + P 2 ) 


(A20) 


Equations (Al 6 ) and (A19) may be compared, with the same object in mind 
as was considered in case I, to show that they are not compatible, and 
therefore that solutions corresponding to values of a unequal to zero 
are not possible for case II. 
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TABLE I.- EXAMPLES OF CASES PLOTTED IN FIGURE l4(a) FOR CASE I 
(J y = CONSTANT), (N 0 = 1000 ) 
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Case 

S 2 

K 2 

7 Mq 2 

7 - 1 ^ 
7 * 

K 2 Bi 

To 

1 

10 2 

10 3 

1 

0.10 

0.39832 

0. 49580x10" 1 




1 

.25 

.15933 

.19832X10- 1 




5 

.10 

.079664 

• 9 9162 X 10 " 2 




5 

•25 

.031865 

• 39664 xio -2 

2 

10 2 

10 

l 

.10 

.39624 

. 49064 x 10” 1 




1 

.25 

.15850 

.19626XIO' 1 




5 

.10 

.079248 

.98129x10“ 2 




5 

.25 

.031699 

.3925IX10" 2 

3 

1.0 

10 3 

1 

.10 

28.288 

2.5007 




1 

.25 

11.315 

1.0003 




5 

.10 

5.6577 

•50014 




5 

.25 

2.2631 

.20006 

4 

1.0 

10 

l 

.10 

21.046 

1.3842 




1 

.25 

8.4186 

.55369 




5 

.10 

4.2093 

.27684 




5 

•25 

1.6837 

.11074 


TABLE II.- EXAMPLES OF CASES PLOTTED IN FIGURE l4(b) FOR CASE II 

(Ey = CONSTANT) 


Case 

§6 

K 2 

No 

7Mo 2 

£-^Pr 

7 

Bn 

T w 

To 

5 

1.0 

10 

1000 

1.0 

.10 

5.3889 

1.4520 





1.0 

•25 

2 .1556 

.58082 




100 

1.0 

.10 

5.3811 

1.4478 





1.0 

.25 

2.1548 

.58038 

6 

1.0 

100 

1000 

1.0 

.10 

8.2955 

3 . 44 o 8 





1.0 

.25 

3.3188 

1.3768 




100 

1.0 

.10 

8 . ion 

3 . 28 l 4 





1.0 

.25 

3.2936 

1.3560 

7 

1.0 

1000 

1000 

1.0 

.10 

9.^137 

4.4309 





1.0 

•25 

3.7726 

1.7791 




100 

1.0 

.10 

7.4200 

2.7528 





1.0 

•25 

3.4657 

1.5014 
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Figure !•- Coordinate system. 
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Uniform applied 
magnetic field, H 0 


(a) Case I (j y = constant; Ey to be determined) . 


Uniform applied 
magnetic field, H 0 



(b) Case II (Ey = constant; Jy to be determined) . 
Figure 2.- Configuration of electromagnetic fields. 
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Figure 3»~ Fully developed flow velocity profiles. 
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(a) Case I. 

Figure 4.- Fully developed flow center-line velocity 
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(b) Case II. 
Figure 4.- Concluded 
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(a) Case I. 


Figure 5*- Temperature in fully developed flow. 
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(b) Case II. 


Figure 5.- Concluded. 
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Figure 6.- Examples of development of flow from uniform profile; case I 
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Figure 6.- Concluded 
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(a) N 0 = 1000, K 2 = 1000, Bjj = 100. 



Figure 7 .- Examples of development of flow from uniform profile; case II. 
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(c) N 0 = 1000, Y? = 1000, B n = 10. 
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(d) N 0 = 1000, K 2 = 100, Bjj = 10. 
Figure 7*- Concluded, 
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(a) Case I. 

Figure 8.- Center-line velocity in accelerating flow relative to inlet velocity. 
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(b) Case II with N 0 = 1000. 


Figure 8.- Concluded 
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(a) Effect of varying K 2 and Bi for case I with N 0 = 1000. 
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(h) Effect of varying No for case I. 

Figure 9*- Center-line velocity in accelerating flow relative to fully 

developed flow velocity. 
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(c) Effect of varying K 2 and Bji for case II with N 0 = 1000. 
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(d) Effect of varying N 0 for case II with Bn = 100. 
Figure 9*- Concluded. 
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(a) Case I with N 0 in the range 250 to 1000. 


Figure 10.- Pipe length required to achieve velocity within 1 percent of 

fully developed flow velocity. 


on ct\ vo 



On vo u> > 


63 


®( 



N 0 

(b) Case II. 


Figure 10.- Concluded 
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(a) Case I. 

Figure 11.- Center- line pressure in accelerating flow. 
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(b) Case II. 


Figure 11.- Concluded 
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(a) Case I. 


(t>) Case II. 


Figure 12.- 


Inlet to fully developed flow pressure difference. 
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(b) Case II (Ey = constant) . 


Figure 13*- Concluded 
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(a) Case I ( Jy = constant) . 
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(t>) Case II (Ey = constant) . 
Figure 14.- Mass flow comparison. 
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